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Abstract: We apply exact WKB methods to the study of the partition function of pure 
M = 2 Cj-deformed gauge theory in four dimensions in the context of the 2 d/ 4 d correspon¬ 
dence. We study the partition function at leading order in 62/61 (i.e. at large central charge) 
and in an expansion in 61. We hnd corrections of the form ~ exp[— periods j expan¬ 

sion. We attribute these to the exchange of the order of summation over gauge instanton 
number and over powers of 61 when passing from the Nekrasov form of the partition function 
to the topological string theory inspired form. We conjecture that such corrections should 
be computable from a worldsheet perspective on the partition function. Our results fol¬ 
low upon the determination of the Stokes graphs associated to the Mathieu equation with 
complex parameters and the application of exact WKB techniques to compute the Mathieu 
characteristic exponent. 
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1 Introduction 

Non-perturbative completions of string perturbation theory have been proposed in various 
backgrounds. One prominent avenue to advance beyond perturbative results is via the 
holographic correspondence between gauge theories and theories of gravity with a negative 
cosmological constant. Another approach (with holographic applications) which has proved 
fruitful is to map string theories or some of their observables to matrix models and to 
compute non-perturbative corrections using techniques such as localization. Nevertheless, 
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the non-perturbative definition of string theory in generic backgronnds remains a wide open 
problem. 

Topological string theory is a promising framework within which to tackle this problem, as 
the natnre of the pertnrbation theory is the same as that of the fnll string theory - giving 
rise to a generically non-convergent genus expansion explicitly linked to computations on 
underlying Riemann surfaces. Any indication as to what type of approximation to a non- 
perturbative theory gives rise to such a perturbation theory could help clarify the structures 
underlying string theory. For a large class of backgrounds, the topological string theory par¬ 
tition function Z^op is fully computable in various series expansions. Lifting these expansions 
to analytic functions would take us a long way towards a non-perturbative understanding of 
the theory. 

The gauge theory limit of certain topological string theories is (ei, e 2 )-deformed Af = 2 
supersymmetric Yang-Mills theories in four dimensions. The genus expansion parameter 
Qs of the topological string is encoded in the ej-parameters via gl = 6162 . The problem of 
determining the instanton partition function Zinst of these theories, which descends from the 
topological string partition function Ztop, has been solved in [1]. The solution is presented 
as a power series in the instanton counting parameter A, with coefficients that are rational 
functions of Ci and 62 . The parameter A reflects the coupling of the gauge theory, which 
maps to certain Kahler parameters in the string theory setting. The disadvantage of the 
resummation of the e^-series is that it occurs at the expense of introducing this new expansion. 
For instance, modular properties of the expansion coefficients of the e^-series are masked in 
the A-series. 

Yang-Mills theories with M = 2 supersymmetries are related to two-dimensional conformal 
field theory via the two-dimensional / four-dimensional correspondence |2]. The powerful 
computational techniques that exist in the framework of two-dimensional conformal field 
theory can thus be put to use to elucidate Cj-deformed gauge theories, and by extension, 
topological string theory. 

The conformal field theory technique at the heart of the analysis in this paper relies on null 
vector decoupling. It permits the computation of conformal blocks, mapped to Zinst under 
the 2d/4d correspondence, via solution of a differential equation. It has been shown [SI IH E] 
that the WKB analysis of this equation, in the limit €2 —)■ 0, reproduces the non-convergent 
ei-expansion of Zinst- Methods exist to enhance WKB results non-perturbatively (see e.g. 
[SI 0 El El HDl El fil DSl El DS] ). The point of departure of this paper is to ask what these 
methods can teach us about the nature of the Cj-expansion of the instanton partition function 
Zinst and eventually the topological string partition function Z^p- Indeed, our analysis will 
yield corrections to this expansion in the form ~ e o. As the instanton partition function 
Zinst in the formulation of [T] is thought to converge [IS], we conclude that these non- 
perturbative corrections arise when the order of summation over powers of the instanton 
counting parameter A and over powers of ei is inverted. At the locus ei = — 62 , we have a 
worldsheet description of the topological string theory, and it gives rise to a (non-convergent) 
expansion of Zinst with this reversed order of summation. Once a worldsheet description of 
the general Cj-deformed theory is formulated (see [m ESI ElEniED for attempts in this 
direction), one should seek to compute the non-perturbative corrections in ei within that 
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framework. Here, we determine them by conformal field theory and exact WKB methods. 
As Zinst in an e^-expansion remains divergent when restricted to the conventional topological 
string locus ei = — 62 , similar non-perturbative corrections in e* should also arise in this more 
standard setting. 

The theory we shall study is pure M = 2 gauge theory in four dimensions, whose instanton 
partition function maps to an irregular conformal block in two-dimensional conformal held 
theory [2^ . The relevant null vector decoupling equation maps to the Mathieu equation with 
complex parameters. We will perform an exact WKB analysis of this equation, determining 
the Stokes regions for various complex values of the parameters, and incorporating Stokes 
phenomena in the computation of the characteristic exponent of its solutions. This proce¬ 
dure introduces corrections of the order exp [—A] in the relation between the characteristic 
exponent, linked to the vacuum expectation value a of the scalar held in gauge theory, and 
a certain complex parameter u of the equation, which is related to the gauge theory parti¬ 
tion function via a Matone-style relation. Inverting this relation maps these corrections to 
non-perturbative corrections to the ei-expansion of the partition function as a function of 
the expectation value a. 


For other studies of non-perturbative ehects in topological string theory, see e.g. [23] . 
where non-perturbative contributions to the spherical partition function of ABJM theory 
are mapped to the topological string amplitude on x P^ together with an additional 
contribution from the rehned topological string amplitude on the same manifold in the 
Nekrasov-Shatashvili limit [ 23 ] . The authors propose this combination as a non-perturbative 
dehnition of the topological string amplitude on this manifold, and conjecture a generaliza¬ 
tion to arbitrary local Calabi-Yau backgrounds. In [ 23 ], the authors consider the topological 
string partition function in Gopakumar-Vafa form, closely related to the Nekrasov form of 
Zjjist lifted to topological string theory, in light of resummation techniques. In I2S1. the 
holomorphic anomaly equations are conjectured to hold non-perturbatively, and invoked to 
conjecture structural properties of a transseries expansion of the topological string free am¬ 
plitude. These ideas are rehned and tested on the example of the topological string on local 
CP^ in [ 27 ], based on computational results for the topological string free energies Fg with g 
up to ~ 100. A similar analysis for the spherical partition function of W = 2 super conformal 
gauge theory and J\f = 2 * gauge theory was performed in [2H]- The paper [ 29 ] studies the 
Mathieu equation in the context of exact WKB and the 2d/4d correspondence, as we shall do 
in the following; unlike the present paper, it restricts to m G M, thus centering the discussion 
around the band structure of the equation that is specihc to this domain. 


The paper is structured as follows. In section [^ we review the relation between the instan¬ 
ton partition function of Cj-deformed pure M = 2 super Yang-Mills theory and irregular 
conformal blocks. We recall the null vector decoupling equation the latter satisfy and map 
it to a standard form of the Mathieu equation. Section is dedicated to the study of this 
equation. Some known facts about the Mathieu equation are summarized in section [3.1 In 
section 3Y, we study the WKB approximation to the solution of the Mathieu equation. We 
go beyond perturbation theory in section |3.3[ We discuss the Stokes graphs associated to 


the Mathieu equation with complex parameters and use exact WKB methods to determine 
properties of analytic solutions to the differential equation. In section we discuss the non- 
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perturbative corrections to the exact periodicity of the irregular block, and its consequences 
for non-perturbative corrections to the e-expansion of the instanton partition function. We 
conclude and list interesting open problems in section Some technical details are relegated 
to appendix while appendices and are dedicated to numerical results. 


2 Pure M = 2 SYM and Conformal Field Theory 

In this section, we review the relation between the gauge theory instanton partition func¬ 
tion of pure M = 2 SYM and the irregular conformal block, as predicted by the 2d/4d 
correspondence [a [22]. 


2.1 The Seiberg-Witten theory of pure M = 2 SYM 


The gauge theory we will be concerned with in this paper is pure M = 2 SU{2) super 
Yang-Mills theory, i.e. the theory of a single vector multiplet with gauge group SU{2). 
This is the original theory solved by Seiberg and Witten [201 using geometric methods. The 
prepotential of the theory is a function of the vacuum expectation value a of the scalar in 
the vector multiplet. The coefficients of all instanton contributions were determined in [1] 
in a form amenable to direct comparison with conformal held theory. 

The Seiberg-Witten curve in |30| was given in the form|^ 

= {x^-l){x-u), (2.1) 


with Seiberg-Witten differential 


A = 


y/x — U 
\/x‘^ — 1 


dx. 


( 2 . 2 ) 


We can choose branch cuts from a; = —ltox = lon the real axis, and from x = u to inhnity, 
not crossing the real axis. A basis of cycles on the Riemann surface (2.1) is then given by 
a curve circling the line connecting x = —1 and a; = 1 in positive orientation without 
intersecting the second branch cut, and the curve circling the points x = 1 and x = u, 
intersecting both branch cuts once and intersecting the hrst curve with intersection number 
1. These choices dehne the A- and the R-cycle of the torus respectively. The corresponding 
Seiberg-Witten periods are 


,(o) _ 


= <p X, = 


A 


(2.3) 


B 


The superscripts indicate that these periods constitute the leading terms in formal power 
series to be introduced shortly. 


Based in part on previous work [211E2], Gaiotto in [22] suggested shifting the emphasis onto 
quadratic differentials when analyzing M = 2 theories. In the case at hand, the quadratic 

mass scale A has here been set to 1. It can easily be introduce by dimensional analysis, with x and u 
carrying mass dimension 2. 
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differential reads 


X — u 


(2.4) 


— 1 


dx ® dx. 


For our upcoming analysis, the variable redefinition 


X = cos q 


(2.5) 


will prove useful, with regard to which (^sw takes the form 


(j>sw = {u — cos q) dq®dq. 


( 2 . 6 ) 


The relation between the variables x and q is one-to-one if we restrict the variable q to the 
range Reg G [—7r,7r], Img > 0 - we will refer to this region as the fundamental domain of 
the g-plane - and identify the half-lines Reg = —vr and Reg = tt, as well as the intervals 
[—TT, 0] and [0, tt] via the map g h-)■ —g. See hgure The branch cut connecting the points 
X = ±1 maps to the intervals [—vr,0] and [0,7r] of the g-plane, undoing the identification of 
these two intervals (as the factor \/l — x"^ = sing in the curve variable y differs by a sign 
between previously paired points). The second branch cut runs from the preimage g^ of u 
in the fundamental domain to positive imaginary infinity, without crossing the imaginary 
axis. We can identify the second sheet of the x-plane with the image of the fundamental 
domain of the g-plane under the map g h->■ —g. Overall, the Riemann surface (2.1) in the g 


coordinate is hence given by a cylinder with branch cuts emanating from g^ and —g^ towards 
the two ends of the cylinder. In this representation of the Riemann surface, the A-cycle on 
the g-plane is represented by a curve running from 7r-l-i|go| to — 7r-l-i|go| without intersecting 
the second branch cut, or simply by a cycle of the cylinder. The R-cycle is represented by 
a curve running from a point on the second branch cut in the fundamental domain of the 
g-plane to its image point in the lower half-plane, crossing both branch cuts once. 

A third choice of variable, underlying the analysis in [221IM] relating pure M = 2 SU (2) 
gauge theory to conformal field theory, is given by 


z = e*" 


(2.7) 


In terms of z, the quadratic differential is (up to an overall factor, and with the scale A 


re-introduced, as it will come in handy in subsection 2.2) given by 


'A^ 2u A^\ 

^SW = I —-V H - I dz®dz. 




z j 


( 2 . 8 ) 


2.2 The irregular conformal block 

The original 2d/4d correspondence |2] conjectured a relation between the instanton partition 
function of super conformal M = 2 gauge theories and Virasoro conformal blocks. This 
relation was extended to asymptotically free gauge theories in [22] by introducing irregular 
conformal blocks. These are defined as the norm of so-called Gaiotto vectors: formal power 
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Figure 1: The q-plane. The 
solid black line is mapped to the 
real axis of the x-plane. The 
interval [—tt, tt] is the double 
cover of the branch cut [—1,1] 
in the x-plane. 


series |A, A^) in the parameter h? with coefficients valued in the Verma module of highest 
weight A satisfying the relations [22] 

Li|A,A2) = A2|A,A2) , L2|A,A2)=0 (2.9) 


order by order in A^. The requirements on the states | A, A^) are chosen such that the ensuing 
expectation value of the energy momentum tensor in the semi-classical limit reproduces the 
quadratic differential (2.8), following the prescription of [2], further elucidated in 


Gaiotto states can be constructed via collision of primary helds |3l]. In this manner, the 
irregular conformal block for pure SU{2) can be obtained by a limiting procedure from the 
conformal block on the sphere with four insertions [31]. Starting from the conformal Ward 
identity satisfied by the product of two primaries acting on the vacuum |0), 


r>(z)TA2(^2)^'Ai(2:i)|0) 


{z - Z2y 


+ ^—d2 + 

Z - Z2 


Ai 

{Z - Zi)2 


+ —^— di) X 
Z — Zi 


'rA,(^2)>tA,(2l)|0), 


( 2 . 10 ) 


where T> indicates the sum over modes of the energy-momentum tensor T with weight larger 
or equal than —1, we can send the insertion zi —)■ 0 to generate a primary state |Ai). We 
moreover use the fact that L_i = 81 + 82 generates a translation on the correlator, as well 
as the parameterizations 


Q = b + b-^ 

Cl = —^ 20^1 

Q. = Oil H" 0^2 1 


( 2 . 11 ) 
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to find, in the second limit Z 2 0 and a and ci fixed, 

2ci{Q — a) a{Q — a) 


r>W'tA,(.-2)'i'A.(zi)|0) = (-:i + 


+ 


Z 2^0 


H— i^cidc^ H—-h-i) X 


( 2 . 12 ) 


To relate this expression to the quadratic differential (2.8) with third order pole, we take 


the further limit A 2 = —cf —?• 0 while Ai = 2ci{Q — a) is kept finite. The state that this 
limiting procedure gives rise to satisfies the Whittaker properties 


L„|A,{W}) =A„|A,{A,}) 
with Ai = A^ and A„ = 0 for n > 1 


(2.13) 


By explicit calculation to a fixed order, it was demonstrated 122 ] that the norm of this state 
coincides with the instanton partition function of Cj-deformed pure M = 2 super Yang-Mills 
theory 

(A, A| A, A) = ^2.14) 

with the conventional 2d/4d identifications 


62 


b = - , p = 


^cft 


61 


27iiy/e^ ’ 


(2.16) 


where the momentum p is related to the conformal weight A by the formula A{ip) = ^+p‘^- 
We will introduce numerous variables closely related to the vacuum expectation value of the 
adjoint scalar in the gauge theory, and have therefore denoted the conformal field theory 
quantity determining the weight of the exchanged state with the subscript ‘eft’. 


2.3 The null vector decoupling equation 

Our point of departure will be the fact that the irregular block after insertion of a degenerate 
operator satisfies a null vector decoupling equation. We denote the irregular block with 
insertion ^h 2 ,i fhe point 2 : as ^/(z). 




(2.16) 


In this definition, we slightly shift the momenta of the incoming and outgoing irregular 


vectors away from 


Q'cft 


27ri^ei€2 


such that 


. / ^eft ^ ^ A f b , 

27ii^€i€2 4 27iiy^ei€2 4 ’ 


^eft 

27riei 


(2.17) 


This condition is necessary for the correlator not to vanish. The derivation of the null vector 
decoupling equation this correlator satisfies was worked out in e.g. [Ml [37]. After a change 
of variables z = e'^ and the rescaling 


'k(z) = 2 ; 


A{ 


— 'Y ' 2'Ki./e-\ eo 


^eft 


62 ^ 4 


+ \)—h 2 . 


W(n;), 


|)-A( 


(2.18) 



















the equation reads 


[-dl + ^ 
£2 


+ —{e^ + e-n + yaA)yH = 0 . 
vr?e2 eie2 4 


To simplify the equation further, we dehne 


‘^cft 


l{w) = ® e 


logA 


Y{w). 


(2.19) 


( 2 . 20 ) 


As the monodromy of ^( 2 ;) that follows from its dehnition (2.16) is accounted for by the 


prefactor in (2.18), the function S(tc) has periodicity 'E{w — 2m) = e o S(tc). It satishes 
the equation 

. A, 


:-5^ + 

£2 


(e-+ e-) +-^a) S(n;) = 0 
6162 4 


( 2 . 21 ) 


This differential equation for S(w) is exact, both in ei and in 62 . 


2.4 The semi-classical limit in the central charge 


To be able to extract the conformal block (2.14) of interest from solutions of the null vector 


decoupling equation (2.21), we will study it to leading order in 62 / 61 . This is a hrst WKB 
approximation. The limit renders the Gaiotto state heavy compared to the light degenerate 
insertion, justifying the factorization ansatz 


l{w) = en' 2 '^*'^^y(t(;|A) , 


( 2 . 22 ) 


with logx(i/;|A) hnite in the limit 62 —)■ 0. The hrst factor is to be identihed to the irregular 
block without the degenerate insertion, while the second factor is associated to the light 


degenerate mode. Plugging the ansatz into equation (2.21), we obtain 

A2 


:-di + 

£2 


(e"' + e "') + —9 a + -9iogA9^)x('w^l^) — 0 . 

6162 4 46162 


(2.23) 


We multiply the equation by 6162 /A^ to hnd 


+ (e'^ + e "') + ^^9iogA9^)x('W^I^) — 0 . 


(2.24) 


The differential equation is exact in 61 . It was already studied (pre-2d/4d correspondence!) 
in [3H] to obtain the prepotential of pure AT = 2 SU{2) gauge theory. Upon the variable 
redehnition 

4 2 1 

w = ^q, ^ = 6 , 

it is mapped to the form 


8A2 


^logA-T U , 


(2.25) 


(-6 9: + cos q)^jJ{q) = u'A{q) . 


(2.26) 


This is the Mathieu equation, which has been studied from various angles in the mathematics 
literature (see e.g. [391110]). 
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3 The Mathieu Equation and Exact WKB 


The technical core of this paper is the exact WKB analysis of the Mathieu equation. This 
method will permit us to compute corrections to the formal WKB solutions to the equation, 
in a sense which we shall make precise in this section. Recall that the formal solutions 
reproduce the non-convergent e-expansion of the instanton partition function of pure M = 2 
gauge theory at large a/e. The correction terms computed here yield information beyond 
this asymptotic expansion. 


Before turning to the WKB analysis in subsection 3.2, we collect some general results re¬ 


garding the Mathieu equation and its solutions in the following subsection, based on the 
exposition in 


3.1 The Mathieu equation 


The Mathieu equation (2.26) is an ordinary differential equation of degree 2. It has a unique 


solution upon specifying two boundary conditions on the solution -0, e.g. the values of 'ip and 
its derivative ip' at a point. The solution is entire in the variable q as well as the parameters 
ue~‘^ and [39]. As the potential term in the Mathieu equation is periodic (of period 27r 
in the variable g), Floquet theory applies. In particular, two independent solutions of the 
equation exist that satisfy 

) («) = , 


(3.1) 


for appropriate z/, with y a 27 r periodic vector valued function of q. Such solutions are called 
Floquet solutions. The factor which reflects the monodromy of the solution ip under q —)■ 
g-|- 27 r, is called a characteristic multiplier of the equation, and z/ is a characteristic exponent. 
Since the characteristic exponent is dehned via the log of the characteristic multiplier, it is 


only dehned modulo 2Z. In equation (3.1), shifting the characteristic exponent z/ by 2Z 
requires rescaling the periodic function y(g). We introduce the parameter z> by writing 
u = 0 + 2n for n G Z and Re (z>) G [—1,1]. The characteristic exponents z) of the two 
independent Floquet solutions of the Mathieu equation add to 0, by the absence of a hrst 


order derivative term in (2.26). We can express the Floquet solutions as linear combinations 


ip = Atpi + Bip 2 of two solutions 'ipi,'ip 2 satisfying the boundary conditions 


V’l ^^2 

^p[ 


( 0 ) = 


1 0 
0 1 


(3.2) 


Requiring that coefficients A and B exist such that the corresponding linear combination sat- 
ishes the Floquet monodromy condition with characteristic exponent u entails the following 
constraint on the characteristic exponent: 


TT, 


COS TTU = 1pl{u, e, —). 


(3.3) 


This equation is called the characteristic equation of the Mathieu equation, or somewhat 
less fortuitously, a quantization condition. For hxed parameters u and e, equation (3.3) 
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determines the characteristic exponent up to a sign and 2Z ambiguity. As we have seen, the 


2Z ambiguity corresponds to a rescaling of the periodic function y;(g) in (3.1) (or a relabeling 
of the Fourier coefficients of 'ijj), whereas the two signs correspond to the two independent 
Floquet solutions. 

In applications, it is often important to know the possible values of the parameter u that are 
consistent with given values of the parameter e and the characteristic exponent v. When the 
Mathieu equation arises as a Schrodinger equation in quantum mechanics, hxing v corre¬ 
sponds to hxing the periodicity condition on the wave function, and determining the values 
of u that permit solutions with this periodicity is tantamount to determining the spectrum 
of the Hamiltonian. In the context of this paper, v corresponds to the exchanged momentum 
of the conformal block S(tc) in (2.20) via 


IVTT = 


^cft 


(3.4) 


By the 2d/4d correspondence, the characteristic exponent hence maps to the scalar vacuum 
expectation value of the adjoint scalar in the vector multiplet of the SU{2) gauge theory. The 
corresponding u determines the partition func tion of the gauge theory via the generalized 
Matone relation mnnES] as it arises in ^2.25^, u = —We will review in subsection 


3.2 that the characteristic exponent at hxed u is approximated by the A-period of the Seiberg- 
Witten differential in a WKB analysis of the Mathieu equation, and study how exact WKB 
methods permit determining corrections to this relation in section 3.3 In the rest of this 


subsection, we will review what can be learned about the relation between v and u from the 


study of the characteristic equation (3.3), following the classic reference 


For a given e and a non-integer z/, a discrete inhnite number of solutions of the equation (3.3) 
for u exists. We can label these as u{u,N,e), with N E h. Notice that by the periodicity 
of the cosine function, we can dehne an integer-valued function N{Mi, M 2 ) on Z x Z such 
that u{u -|- 2Mi, M 2 , e) = m(z/, N{Mi, M 2 ), e) for Mi, M 2 G Z. One can moreover show that 
the choice N{Mi, M 2 ) = 2(Mi -|- M 2 ) is possible, allowing us to combine the variables u and 
N and express the solution u as a function m(z/ -|- 2N,e). No generality is lost by calling 
the hrst variable u. Recall that by Floquet theory, we furthermore have the parity property 
m(z/, e) = u{—u, e). 

At integer z/ = n G Z — {0}, two discrete inhnite families of solutions to the characteristic 


equation (3.3) exist, labelled as u+{n,e) and U-{n,e) (these solutions, rescaled by the factor 
4e“^, are usually denoted a„ and bn). The function m(z/, e) for z/ G M is discontinuous at 
z/ G Z, 

lim u{h',e) = u±{n,e) , (3.5) 

with z/ —)■ rz,^ indicating that rz G Z is approached from above/below respectively. These 
discontinuities along the real axis give rise to what is referred to as the band structure of 
the spectrum of the Mathieu equation, see hgure 

Convergent and asymptotic series in e for u are known, at large and small e respectively. It 
can be uniquely characterized by its convergent power series expansion around e = 00 , 




4 4(z/2 — 1) 




(3.6) 
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Figure 2: The variable u is discontinuous as a function of real v. On the left, the behavior 
% ~ is visible, as well as the narrowing of the bands as e —>■ cx). On the right, the limit 
M —)■ —1 as e —)■ 0 is apparent. 


The coefficients of the series expansion can be found by equating negative powers of e in the 
following continued-fraction equation. With the notation x = 4ue~‘^, 


,-2 


,-2 


,-2 


,-2 


_ 2 _ _ _ _ _ _ _ _ _ _ _ fo 

^ ^ X — {u + 2y—X — {u + 4^— X — (z/— 2)2—X — (z/ — 4)2— 

The series that results from this procedure has hnite convergence radius when the parameter 
u is not an integer. For integer z/, convergent power series expansions around e = oo exist 
as well, with leading term 

1 1 

(3.8) 


J + 0(i), 


but the corresponding coefficients are given by recourse to different continued-fraction equa¬ 
tions. E.g., for M_|_ and even n. 


,-2 


,-2 


,-2 


X — (2n)^ — 


2 e 


-2 


X — {2n — 2)2— X — (2n — 4)2 — 


x-22- 

^-2 


X 


,-2 


(2rz,-I-2)2 — X—(2n-I-4)2 — X— ' 

Notice that this continued-fraction equation is different in nature from the one appearing in 


equation (3.7), as the variable n also appears as a summation index. Hence, the equation 
yields the coefficients of the series expansion only once n has been assigned a value. For 
large n, the leading coefficients in the series expansion can nevertheless be expressed as 
rational functions of the parameter rz,, and they coincide with the coefficients for v non¬ 
integer, evaluated at z/ = n. This agreement breaks down for higher terms in the series. 


rendering (3.8) convergent, whereas (3.6) diverges for integer v 


For real and small, it is known that the values of the two solutions 'U+(n, e) and M_(n-Tl, e) 
approach each other and have the asymptotic expansion, for positive. 


M+(rz,e) ~ M_(n + l,e)-1 + + l)e^ + ... , 


(3.10) 
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where s = 2n + 1, n = 0,1, 2,.... A similar expansion can be derived for small and 
negative, with leading term +1. One can show that the diffe rence between u+{n,e) and 

U-{n + l,e) scales like exp The asymptotic expansion (3.10) is hence also valid for 

m(z/, e), with u real and in the interval n < u < n + 1. 

Finally, one can ask about the analytic properties of u{v, e) as a function of v. According to 


, the solution 'ipi on the right hand side of the characteristic equation (3.3) is analytic as 


a function of v and e away from possible branch cuts. By the implicit function theorem, the 
parameter u will be analytic as a function of v and e away from these branch cuts and from 
zeros of the w-derivative of F(m, z/) = "01 (m, e, |) — costtz/. By writing 


f) 

du 


dv 

—TT— sinTTZ/ , 
ou 


(3.11) 


the discontinuities of u for real z/ become visible as the zeros of sinvrz/. Additional analytic 
complications will arise at zeros of the factor dui'- Using Mathematica, we found some 
evidence for the existence of one stationary point UQ{n) of z/(m) for each pair u±{n), rz, G Z, 
with u_{n) < Uo{n) < M+(n). 


3.2 WKB analysis of the Mathieu equation 

In this section, we will review the WKB approximation to the Mathieu equation. This is 
standard material. For future convenience, we use the conventions in [T3]. For further results 
on the WKB analysis of this equation, see Hami]. 

The starting point of the general theory is a second order differential equation of the 
form 

[62aJ-Q(g,6)]vI/(g) = 0, (3.12) 

on a Riemann surface S, depending on a potential Q{q,e) = which is a poly¬ 

nomial in e with coefficients Qk{q) that are meromorphic on S, satisfying conditions outlined 
in [13]. We will immediately specialize to the Mathieu equation in the form 

— (cosg — M)]T(g) = 0 , (3.13) 

with S chosen as the cylinder —vr < Im g < tt compactified to a sphere by adding points at 
±ioo. The parameters e and u can be complex. Hence, 

Qo = cos q — u (3.14) 

and, since we pick the parameter u to be e-independent, we have Qn = 0 for n > 0. 

The WKB ansatz for the solution of the differential equation is 

'ip{q,e) = exp{j S dq), (3.15) 
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with S expanded as a formal power series in e, 


S=^S_i + So + eS, + ... . (3.16) 

By plugging this ansatz into the differential equation, one immediately derives the following 
recursion relation for the coefficients S'„: 


2S'_iS'n+i + SniSn2 + 


ni-\-n2=n 


Si, 

dSn 

dq 


Qo , 

Q n+2 


n > —1. 


(3.17) 

(3.18) 


Equation (3.17) has two solutions, = ±y/Qo. We note that up to normalization, the 


Seiberg-Witten differential of pure Af = 2 SU{2) gauge theory thus arises in the WKB 
analysis of the Mathieu equation as A = S^idq. In this matching, the variable u defined 


via (2.25) coincides with the conventional variable parametrizing the Seiberg-Witten w-plane 


introduced in subsection 2.1, a manifestation of a generalized Matone relation. 


Depending on the choice of the sign of S^,, we obtain two solutions to the recursion relations, 
which we label S^. The first few series coefficients, for u chosen to be e-independent, are 
given by 

= ±a/cos q — u , 


S^ 

Sf 




sf 


--d\ogS-i dq = - -, 

2 4 cos q — u 

^ cos 2q + 8m cos g — 9 

64(cosg — m)^/2 

sin(g) (20Mcos(g) + cos(2g) + 8 m^ — 29) 


128(m — cos(g))^ 


= T 


16384(cos(g) 


(iQ {32u^ - 265) M cosg + 20 - 173) cos2g + 


912m cos 3g + 25 cos 4g - 1344m^ + 5355^ . 


Introducing 


one can show that 


Sodd = ^( 5 + - ^-) , = ^(^+ + ^-) , 

^ _ 1 d log Sodd 

^even — ^ 7 

2 dq 


The formal WKB solution (3.15) thus takes the form 

1 


^/'±(?,e) = 


\/ Soddiq^ *^) 


exp ± / Soddiq,e) dq , 


(3.19) 

(3.20) 

(3.21) 
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where this formal expression is to be interpreted as an analytic function in q with branch 
cuts multiplying a formal power series as follows: 


= exp ( ±- 




(3.22) 


k=0 


To £x the normalization of (3.21), we need to specify the starting point of the integral. In 
this paper, we will choose this starting point to coincide with the zeros of Qo, called turning 
points. Care is required in dehning the ensuing integral [H], as the coefficients Sn, n > 0, 
have poles at the turning points. 

Our analysis will require comparing WKB solutions normalized with regard to different 
turning points. These are related by exponentials of periods of Sodd, for which we introduce 
the notation 


a = e 


Sodd dq , Qif) 6 / Sodd dq ■ 


(3.23) 


’ B 


The integrals are to be understood order by order in e. To leading order, the integrals a 
and an equal the Seiberg-Witten periods and which can be expressed in terms of 
hypergeometric functions. 


i(o) = , 

\ 2 2 u + 1J 


.(0) _ 


TT 


''D 




V2 


1 1 


2 2 


(3.24) 

(3.25) 


Instead of evaluating the integrals over the coefficients Sn, n > 0 directly, one can dehne 
differential operators D 2 n of order 2n that map S_i to S 2 n-i up to total derivative terms 
[ami]. The hrst few of these are HD 


D2 

Da 

De 

Ds 


— { 2 udl + du), 

1/28 2o4 8 03 

^^ 45 “ ^n+ ^udn + -dj , 

1 ,124 3 ^g 158 2^5 153 ^4 41 3 . 


(3.26) 


1 / 127 13 3^7 

^^23 X 4725“ 


517 2 ofi 9539 

_7/^(9° -I_ 

24 X 63 “ 23 X 945 


05 , 15229 ^ 4 , 


Derivatives of the hypergeometric functions can in turn be rewritten as hypergeometric 
functions. By acting with the D 2 n on and one hence obtains the coefficients of 
the formal power series (3.23) again in terms of hypergeometric functions of the modulus 

u. 


Note that expressing the higher order corrections to the A- and i?-period via derivative op¬ 
erators acting on the leading contribution establishes that the monodromy matrix governing 
the behavior of and upon circling singularities in the u-plane is not corrected at 
any order in e: schematically, D logu) = logu -|- logu, and D logu does not 
exhibit monodromy. 
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The two formal series allow us to eliminate u, and to solve for instance for J^(a, e), which 
agrees with the perturbative expansion of the Nekrasov instanton partition function. 

We performed successful numerical checks on the perturbative WKB approximation in a 
range of parameters where the Stokes phenomena we shall discuss in the next subsection are 
numerically negligible. We give an example of such a check in appendix |B.1[ 


3.3 Beyond perturbation theory 


We have seen that the monodromy of the solution to the null vector decoupling equation 


(2.24) plays an important role in our analysis: it corresponds to the exchanged momentum 


of the conformal block or, equivalently, to the scalar vacuum expectation value in the gauge 
theory. We also know by Floquet theory that a basis of exact solutions to the Mathieu 
equation (3.13), called Floquet solutions, exists with monodromy behavior ^p± —>■ 

As we have reviewed above, the WKB ansatz gives rise to two independent formal solutions 
of the differential equation as a power series in e which formally diagonalize the monodromy 
matrix; they hence approximate in a sense we shall discuss presently the Floquet solutions 
of the Mathieu equation. Their monodromies under g —)■ g + 2ni to leading order in e are 
given by ± the A-period integral of S-i, which can be identihed with the Seiberg-Witten 
differential A. This period hence provides an approximation to the characteristic exponent 
u via a ~ zTrez/. In this section, we will see how to incorporate exp[— 1/e] corrections in this 
analysis. We introduce the notation for the A-period incorporating these corrections, such 


that = mev. In terms of the quantities introduced in section 2.2 


t/ei- 


The formal power series 'ip± obtained from the WKB ansatz in the form (3.22) generi- 


cally do not converge. Instead, they provide asymptotic expansions to actual solutions 
to the given differential equation: for e contained in a sector arge G ( 6 *i, 02 ) (we will 
have much more to say about this range in the following), solutions Tj- exist of the form 

= exp \/Qo{z) dz^ e^/^/±(g,e), with f± analytic in g and in the given sector for 

e, such that for any N eN and p > 0, a constant C > 0 exists with 


N-l 


\-N 


l/±(g,e) - X] <C Ve : |e| < p,arge G ( 01 , 62 ). 


(3.27) 


k=0 


It is easy to see from the dehnition that a function with an asymptotic expansion in terms 
of formal power series has precisely one such expansion. 

Borel resummation is a technique, given a formal power series in e, to construct a function 
analytic in a sector of the e-plane which has the formal series as its asymptotic expansion. 
The Borel sum is constructed in two steps. The Borel transform 'ipB of a formal series "0 is 
dehned as 

00 00 n —1 

^(e) = X]0fce^ ^ ijBiy) = • (3.28) 

fe =0 k=l ^ '' 

If 05 converges around p = 0 and can be analytically continued along a half-line connecting 
0 to inhnity at an angle —9 to the positive real j/-axis, we can dehne the Laplace integral of 
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V’B in direction 9 as 


Sg['llj]{e) = -00 + 


(3.29) 




'le 


If this integral exists, Se[ij\ provides the sought after analytic function; it is called the Borel 
sum of '0 in direction 9, and 'ip is called Borel summable. Notice that at given e, the existence 
of the integral generically constrains 9 to lie within the sector 


9 e 



— arge 


71 

, 2 ~ ^ 


)• 


(3.30) 


The WKB analysis of the Mathieu equation gives rise to formal WKB series that are as- 
sumec^ to be Borel summable away from a discrete infinite set of angles 9. These angles 
partition the y-plane into sectors, half of which, given an e and in accord with (3.30), de¬ 
termine Borel sums of pj, possibly differing by exponentially suppressed terms amongst each 
other. This ambiguity or integration path dependence of Borel resummation gives rise to 
the so-called Stokes phenomenon. 


When the asymptotic series being resummed is the WKB solution to a differential equation, 
its coefficients depend on a parameter q, and Borel resummation under favorable circum¬ 
stances leads to analytic solutions of the differential equation. For an ordinary differential 
equation of second order, this space is two-dimensional. The Stokes phenomena thus cor¬ 
responds to assigning a different linear combination of a given basis of analytic solutions of 
the differential equation to WKB solutions via the process of Borel resummation. 


When considering g-dependent coefficients, two types of singularities, mobile and fixed, can 
appear in the Borel plane; mobile singularities are those whose position depends on q. Away 
from isolated points in parameter space, all singularities that appear in the Borel plane are 
mobile. Our analysis will hence essentially focus on this case. Keeping the integration path 
of the Laplace transform in the Borel plane fixed, the Stokes phenomenon in this context 
manifests itself by discontinuities in the Borel resummation when q crosses certain lines, 
called Stokes lines, on E. The Stokes lines divide S into domains called Stokes regions. By 
the foregoing, the two solutions of the Mathieu equation that we obtain by Borel resummation 
of the formal WKB solutions depend on the Stokes region: the analytic continuation of the 
Borel resummation into a different Stokes region will equal a linear combination of the Borel 
resummed solutions native to that region. In other words, Borel resummation and analytic 
continuation in q do not commute (before Borel resummation, analytic continuation is to be 
understood term by term). 

Returning now to the question of determining the characteristic exponents of the Mathieu 
equation, we see that two phenomena need to be taken into account when passing from WKB 
to exact results. Firstly, the period of Sodd, which naively coincides with the characteristic 
exponent, must be Borel resummed. The resulting analytic function of e, called a Voros 
multiplierwill depend on the integration path chosen for the Laplace transform. Secondly, 

^The mathematical literature on this subject is uncharacteristically beset by assumptions and deferred 
proofs. But see [45] for a proof of Borel summability in the case of a particular polynomial potential, and 
the forthcoming work [46] for the general polynomial case. 

^In the literature, this term is sometimes also used to indicate the formal period before Borel resummation. 
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to determine the monodromy matrix of a pair of solutions requires analytically continuing 
them between different Stokes regions. As a consequence, we will see that the Floquet 
solutions to the differential equation do not coincide with the Borel resummation of the 
WKB solutions. Both manifestations of the Stokes phenomena must be taken into account 
to determine the characteristic exponents of the Mathieu equation. We will show that 
apparent ambiguities due to the choice of integration path cancel out in the process. 


3.3.1 The Stokes graphs 


Studying the Borel resummation behavior of the WKB solutions requires determining the 


Stokes graphs of the Mathieu equation (3.13) for a given choice of the parameters u and e. 


These are entirely determined by the leading contribution Qq to the potential specifying the 


differential equation (3.12). More precisely, upon a change of variables q ^ q and a rescaling 


of the solution d' to absorb the ensuing hrst order derivative, (3.12) remains form invariant 
upon replacing Qo(z) by 

Qo(S,e) = Qo(z(z) f , (3.31) 


(3.32) 


dz ^ 

and shifting the higher order coefficients Qn- The invariant quantity is hence 

(j) = Qodz^\ 


It is thus natural to interpret Qq as the coefficient of a section of the line bundle 
with K the canonical line bundle on S. Such sections are called quadratic differentials. We 
encountered them in the context of Seiberg-Witten theory in section 2.1 As we remarked 


above, the square root of the quadratic differential, A = S-i dq, coincides with the Seiberg- 


Witten 1-form of the underlying gauge theory; consequently, the quadratic differentials (3.32) 


and (2.6) coincide (up to irrelevant normalization). 


To render A single-valued, we introduce the double cover S of the Riemann surface E, 
branched at the simple zeros and poles (if present) of 0. Note that since we have not 
restricted Reg > 0 in defining S, S does not coincide with the conventional Seiberg-Witten 
curve. In particular, the branch cuts in figure in the upper and lower half-plane are not 
identified on S, and a path connecting them does not yield a cycle. Since the two branch 
cuts can be chosen to be mapped into each other under g —)■ —g, and the formal power series 


Sodd which determines the formal WKB series (3.21) is even in g, we can essentially ignore 


this subtlety until the end of subsection |3.3.2[ In an abuse of terminology, we will continue 
to refer to a path connecting the two turning points as the R-cycle or as homologous to the 
R-cycle. 


We define trajectories of the quadratic differential 0, called WKB curves by [32], by the 
condition that A have constant phase along them. In other words, tangent vectors dt to 
trajectories satisfy 


A ■ a* e M+ . 


(3.33) 
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This translates into the integral condition 


Im / A = const. 


(3.34) 


Distinguished points on the Riemann surface S are given by the zeros of 0. These are called 
turning points of the differential equation. Depending on the order of the zero, we distinguish 
between simple, double, or higher order turning points. A Stokes line is a trajectory that 
ends on a turning point go, hence satishes the equation 


lme^^rX = 0. (3.35) 

Jqo 

The graph formed by these Stokes lines is called the Stokes graph in the direction 9. As 
noted above, the relevance of Stokes lines stems from their relation to the position of the 
poles of the Borel transform of the formal WKB solutions: when q lies on a Stokes line 
associated to an angle 0, the Borel transform exhibits a pole on the line 


For the analysis of transition behavior between Stokes regions, it is important to endow 
Stokes lines with an orientation. We will define the real part of e®® p A to increase in the 
positive direction along a trajectory. This convention implies that a Stokes lines is oriented 
away from the turning point go if Re ^ ^ along it. 


By equation (3.35), we see that the pattern of Stokes lines is determined by 6 * G [0, tt); under 
6 ^ 6 + IT, the pattern remains invariant, but the orientation of each Stokes line flips. 


A local analysis establishes that the number of Stokes lines emanating from a turning point 
is determined by its order, as follows: an order n turning point leads to local behavior 
of the integral (3.34) and therefore has n + 2 Stokes lines emanating from it, with angle 
27r/(n + 2 ) between two neighboring lines. Two types of Stokes lines will be relevant for our 
analysis: simple (or single or separating) Stokes lines, which run between a turning point and 
a pole of 0, and double Stokes lines (or Stokes saddles) that run between two turning points. 
Simple Stokes lines that are oriented away from turning points will be called dominant, those 
oriented towards turning points recessive. The Stokes lines emanating from a turning point 
are alternately dominant and recessive, except upon crossing a branch cut. 


The Stokes graphs of the Mathieu differential equation have been studied in the Gaiotto- 
Witten variables in |32]. We will study them on the g-cylinder more directly related to 
the traditional form of the Mathieu equation. For the computation of non-perturbative 
corrections to the approximation to the characteristic exponent given by the A-period of the 
Seiberg-Witten differential, we will need to study the Stokes graphs for any complex value 
of the parameter pair (m, e). 

In our variables, the quadratic differential 0 is given by 


0 = (cos q — u) dq® dq . 


(3.36) 


Its turning points on the cylinder he at 

Qup/down £ COS 


(3.37) 
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with the subscript indicating the g-halfplane on which the respective preimage of u lies. The 
differential 0 has an essential singularity at the infinity of the complex plane. To analyze 
its behavior restricted to S, it is convenient to revert to the variable z = cosg, to discover 
a pair of cubic poles, one at each point at inhnity on E. The behavior at such poles is that 
trajectories which approach sufficiently are attracted to the pole along a tangent line (more 
generally, the number of tangent lines equals 2 less than the order of the pole) [TF] . 


Critical Stokes graphs Fixing a pair (M,e), studying the occurrence of double Stokes 
lines at special angles 9, called critical, is important for several reasons. We can imagine 
the y-plan^ being divided into sectors via the critical angles. The global topology of Stokes 
graphs is constant within each sector, and transitions upon crossing into a neighboring sector 
in a simple manner. Pairs of sectors related by reflection through the origin differ only in the 
orientation of all Stokes lines. Voros multipliers exhibit jumping behavior upon transition 
between sectors. 

We can make qualitative statements about the occurrence of double Stokes lines by recourse 
to an observation of |1SII32] matching BPS states in the spectrum of the Af = 2 gauge theory 
associated to a given Seiberg-Witten curve with the occurrence of such lines. The charges 
of the BPS state correspond to the homology class associated to the line (recall that the 
double Stokes lines connect turning points on different sheets of the Riemann surface, hence 
determine closed curves on it). From our knowledge of the BPS spectrum of pure M = 2 
gauge theory [301 HO], we are thus led to distinguish two regions on the n-plane, separated 
by the curve of marginal stability, which runs through the points u = ±1. Inside the curve, 
the BPS spectrum consists of the monopole and the dyon, of charge ±(0,1) and ±(1, ±1) 
respectively (the relative sign between the electric and magnetic charge of the dyon is not 
monodromy invariant). Hence, at fixed e and any value of u within this region, two values of 
the angle 9 in the interval [0, tt) should give rise to Stokes graphs exhibiting a double Stokes 
line. One of these will connect the turning points and qdown directly, the other will 
wrap around the cylinder once before connecting the turning points. These angles shifted 
by TT correspond to the respective antiparticles. Outside the curve of marginal stability, the 
spectrum consists of inhnitely many BPS particles, of charge ±(n, 1) for n G Z (in addition 
to the vector bosons at charge ±(1,0), which also give rise to distinct Stokes patterns, see 
[IHl [32]) • At a given value of u in this region, inhnitely many values of 9 in the interval 
[0, tt) hence give rise to double Stokes lines, one for each wrapping number n G Z around 
the cylinder. 

We now turn to a more systematic study of the critical graphs. 

Critical Stokes graphs at u = ±1 A good starting point for the systematic study of the 
Stokes graphs of the Mathieu equation is at the values of u at which the quadratic differential 
exhibits double turning points, i.e. at the monopole and dyon points u = ±1. According 
to our general analysis above, n ± 2 = 4 Stokes lines emanate from each such point. For 
concreteness, let us consider the point u = —1. A critical Stokes graph for this choice of u 

® Recall that y is the Borel dual variable to e. 
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Figure 3: Stokes graphs: the graph in the first row is for u on the curve of marginal stability at 
M = — 1, the second row depicts the two critical graphs inside the curve of marginal stability, 
and the third row a selection of critical graphs outside the curve of marginal stability. 
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occurs at 6 * = 0, as A ■ G M along the line connecting the turning points at g = —tt and 
g = TT along the real axis. This double Stokes line corresponds to the monopole. The dyon 
at M = — 1 is massless, the corresponding double Stokes line has zero length. It arises at 
0 = I. The corresponding two graphs are depicted in the hrst row of hgure 

Critical Stokes graphs: inside the curve of marginal stability Moving away from 
the singular points u = ±1 into the strong coupling region inside the curve of marginal 
stability, each double turning point splits into two single turning points. These have n+2 = 3 
Stokes lines emanating from them. Keeping u real, the two critical Stokes graphs still lie at 
6 * = 0 and 6 = |, as depicted in the second row of hgure The dyon acquires a mass, as 
the double Stokes line corresponding to 0 = | now has hnite length; along it, A is purely 
imaginary. Giving u an imaginary part moves the turning points off of the real axis, while 
maintaining the topology of the diagrams. In particular, the two simple Stokes lines at each 
turning point run off to imaginary inhnity in opposite half-planes. 


Critical Stokes graphs: outside the curve of margiual stability Starting from 
M = — 1 , moving into the weak coupling region outside the curve of marginal stability while 
keeping u real gives rise to the Stokes graphs depicted in the third row of hgure The double 
Stokes lines in the leftmost diagram correspond to the n —?• cx) limit of the BPS particles 
of charge na + qd- Within a small interval around 6 = 0, an inhnite number of double 
Stokes lines arise, which wrap, as 6 approches 0, an increasing number n of times around the 
cylinder before connecting the two turning points, corresponding to central charge na -|- a^i. 
The second graph in the last row of hgure corresponds to the value n = 2. Moving away 
from 6 ^ = 0, n decreases, till it reaches 0 at 6 = |, as depicted in the hnal graph in the 
third row. Further increasing 6 yields double Stokes lines that wrap the cylinder in the 
opposite direction. Giving u an imaginary part shifts the turning points away from the line 
Reg = ± 7 r, while maintaining the topology of the diagram. In particular, aside from the 
two critical graphs corresponding to n —)■ oo with two double Stokes lines attached to each 
turning point, all other critical graphs exhibit two simple Stokes lines at each turning point 
moving oh to imaginary inhnity in the same half-plane. 


In accord with the BPS analysis above, the ?/-plane is hence split into four sectors inside the 
curve of marginal stability. Outside this curve, it is split into inhnitely many sectors that 
accumulate at 6 * = 0 and 6 = tt. We will introduce a convenient indexation of these sectors 
below. 


Generic Stokes graphs As we move oh of a critical value of 6, the double Stokes line £q 
splits into two simple Stokes lines. These can avoid each other in two topologically distinct 
manners: upon decreasing 6 away from a critical value, the lines swerve to the left of io as 
seen from the turning point from which they emerge, upon increasing 6, they swerve to the 
right. This behavior is visible in hgure for u lying outside the curve of marginal stability, 
and again in hgure [^for u inside the curve. The Stokes graphs on the two sides of the critical 
0 -value are said to be related via a hip. 
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Figure 4: Two Stokes graphs related by a flip. If 9q denotes the critical angle, and 5 is a 
small positive constant, the graph to the left is representative of the topology at 6q — 6, the 
graph to the right of that at 9q + 6. 

The generic Stokes graph hence has three Stokes lines emerging from each turning point 
and running off towards inhnity, two towards inhnity in the upper half-plane of the g-plane 
and one towards inhnity in the lower half-plane, or vice versa. Inside the curve of marginal 
stability, these two cases are interchanged by a hip. Outside this curve, the turning point 
in the upper half-plane always exhibits two simple Stokes lines running oh to inhnity in the 
positive half-plane. By the symmetry q —)■ —q of the equation, the behavior at the alternative 
turning point is obtained by rehection through the origin. 


3.3.2 Computing non-perturbative corrections to the characteristic exponent 


Having determined the Stokes graphs of the Mathieu equation, we are now in a position 
to compute contributions to the monodromy of its solutions that are not visible via formal 
WKB analysis. 


As explained in the introduction to this section, the Borel resummations of the two formal 
WKB solutions ^)) at hxed e, yield a diherent basis of solutions of the Mathieu equation 
depending on the Stokes region in which the argument q lies. The solutions obtained 
upon Borel resummation with q in Stokes region (f) can be analytically continued into a 
neighboring Stokes region (j), yielding a basis of solutions also here. A second such basis 
can be obtained directly by Borel resummation of 'ip±{q, e) with q chosen in region (j). The 
matrix relating these two bases of solutions is referred to as the connection matrix 

from region (i) to region (j). It depends on the choice of normalization of the WKB solutions. 
In this paper, we will, as pointed out in section 3.2, normalize the WKB solutions 'ijj± at 
the turning point from which the Stokes line emanates that we wish to cross. Denoting the 
turning points as qk, we will use the notation {qk)± to indicate the Borel resummations in 
Stokes region (i) of the WKB solutions normalized at turning point qk, with q^ lying on a 
boundary of the Stokes region (i). 
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In this notation, 


((rf’) ((«)i) 

where labels a point in Stokes region (j) neighboring Stokes region (i). With the given 
choice of normalization, the matrices S'(,)^q) take a simple form ([8], theorem 2.25 of [H]): 
upon analytically continuing across a dominant Stokes line counterclockwise with regard to 
the turning point, it is given by [H] 


c -MM 

^dom I 0 1 / ' 


(3.39) 


whereas analytic continuation counterclockwise across a recessive Stokes line, requires the 
connection matrix na 

J ) . (3.40) 

As a consistency check, note that a full revolution around a turning point yields the unit 
matrix as connection matrix: 


SrerSi 


rec^branch^rec>^dom 


Sre.r.Sdnm. — 



Sdom S branch Sdom Srec 



(3.41) 


with 


S, 


branch — ^ 



(3.42) 


the transition matrix upon crossing a branch cut. The factor —i is due to the square root 


in the denominator of (3.21). 


From the form of the WKB solution (3.15), we can read off the following relation between 
two solutions (q'i)±^ and {q 2 )±-i when qi and q 2 lie on the boundary of the same Stokes 

(i)-. 


region 


= 


( exp 

r<i2 

± / Sodd 

\ 

- 


{<12 )± 


(i) 


(3.43) 


The prefactor of {q 2 )^± is the Voros multiplier associated to the cycle represented by the line 
connecting the two turning points qi and q 2 '- the notation (■),, denotes the Borel resummation 
of the formal power series in parentheses. The subscript s indicates the sector in the i/-plane 
in which this resummation is performed: as we will review below, Voros multipliers are 
locally constant functions of 0 in a given sector; they can jump as 9 crosses a critical angle. 
We have suppressed the sector dependence elsewhere in the notation. 


The A-monodromy outside the curve of marginal stability Computing the A- 
monodromy for any value of the parameter pair (m, e) with u outside the curve of marginal 
stability, requires crossing at least two Stokes lines. To standardize our calculations, we will 
always start off at a point g in a Stokes region, henceforth Stokes region (1), chosen such that 
the first Stokes line to cross in analytically continuing the solution along the negative A-cycle 
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u=-2 + 2f , 0=0 



Figure 5: The Stokes regions (1), (2), (3) and turning points qi, g 2 , Qs as introduced in the 
text. 


(recall from section 2.1 that the A-cycle runs from tt to — tt in the g-plane) is a Stokes line 
connected to a turning point in the upper half-plane, qi. The analytic continuation across 


this line into Stokes region (2) will involve one of the two connection matrices (3.39) or (3.40), 
depending on whether the Stokes line is dominant or recessive. The next Stokes line to cross 
is a Stokes line emanating from a turning point, g 2 , in the lower half-plane, and has opposite 
orientation. We thus arrive in Stokes region (3), which is identihed with Stokes region (1) 
by the periodicity of the problem. Choosing our branch cuts, as we shall do throughout, to 
connect the turning points in the upper/lower half plane to imaginary inhnity in the same 
half plane (see £gure[^, no branch cut is crossed along this path of analytic continuation. 
The terminology introduced here is exemplihed in hgure 

The terminology introduced in the previous paragraph allows us to introduce the indexing of 
sectors of the |/-plane promised above: given qi G [— 7 r-l- 27 r/cj, 7 r-|- 27 r/cj] for ki ^2 ^ we dehne 
n = fci — ^2, such that the cycle connecting qi and q 2 is homologous to ±{nA + B), with the 
orientation of the cycle chosen to coincide with that of the double Stokes line connecting the 
two turning points at the appropriate boundary of the ^-sector. The sign is correlated with 
whether the first Stokes line crossed is dominant (-I-) or recessive (—). At a given choice 
of (u, e), the sectors of the j/-plane are hence uniquely indexed by (n, ±). If an angle 9 lies 
within the sector [n, ±), then 6 + 7i yields the same Stokes graph with the orientation of all 
Stokes lines reversed, hence lies in the sector (n, =p). In accord with (3.30), exactly one of 


these choices of 9 is compatible with the phase of the parameter e, which we have kept hxed 
throughout the argument. Flipping the sign of e necessitates making the alternative choice. 
In this sense, flipping the sign of e results in the map {n, ±) i—)■ (n, q^). 


e !-)■ -e 


(n,±) ^ (n,T) 


(3.44) 
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Introducing the notation Ni^j for the matrix mapping a basis of solutions obtained via 
Borel resummation of formal WKB functions normalized at turning point g* to the basis 
obtained from WKB functions normalized at turning point the analytic continuation we 
have described is encapsulated in the following equation: 


(gi + 27r) \ \ - AT c 


■ 5 !- 


>( 3 ) 




(b- 



(g + 27r). (3.45) 


By the periodicity of the differential equation, 


(gi + 27r)^^^ 
(gi + 27r)L^^ 


(g + 27r) = 



(3.46) 


thus permitting us to identify the product of matrices in (3.45) as the monodromy matrix 


Ma (recalling again that the A-cycle runs from tt to —tt). With the notation introduced in 


subsection 3.1 and by Floquet theory, 


Tr MA = e^ + e 


(3.47) 


With the conventions introduced above, the appropriate normalization matrices Ni^j are 
given by 


N- 


l-s-2 — 



(3.48) 


and 


aj^-\-{n+l) a 


N, 


2^-3 


0 

ajj + in+1) a 


(3.49) 


The formal periods a and were introduced in (3.23), and the subscript s denotes the 


sector in which the Borel resummation is to be performed. 


Substituting all matrices into equation (3.45), we obtain the following trace of the mon¬ 
odromy matrix: 

Tr Ma = f2 cosh - + eT7Wi+2^)+2'^o)') (3.50) 

V e / {n,±) 

At a given choice of the parameter pair {u, e), a unique answer for Tr Ma exists - in particular, 
this answer cannot depend on the choice of integration direction 6 of the Laplace transform, 
i.e. the choice of sector (n, ±). We hence need to explain the two apparent 6'-dependencies of 
this result: that we have a priori arrived at two different expressions, depending on whether 
the first Stokes line crossed is dominant or recessive, i.e. whether 6 lies in the sector (n, -|-) 
or (n, —), and the n-dependence of these expressions. The key to resolving the first apparent 


ambiguity lies in the relation (3.30): the phase of e determines whether the sector (n, -|-) or 


(n, —) is appropriate. Moreover, by (3.44), flipping the sign of e leaves the result invariant. 
To address the apparent n dependence of the result, we will need to discuss the jumping 


behavior for Voros multipliers. We will do this in subsection 3.3.3 below 
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Figure 6: Inside the curve of marginal stability. 


The A-monodromy inside the curve of marginal stability Representative graphs for 
the two sectors (up to orientation inversion) that arise inside the curve of marginal stability 
are depicted in hgure The graphs in the hrst row arise at Imu > 0, those in the second 
at Imu < 0. The topology of the graphs [^/|^ can be distinguished from that of the graphs 
[^/|^ by whether the turning point in the upper half-plane is connected to imaginary inhnity 

in the upper half-plane via one or two Stokes lines. We introduce the notation (-1-) and 

(-1- H—) to distinguish the corresponding 6*-sectors. The graphs that arise in the sector 
(-1- H—) for Imu > 0, Imu < 0 match the graphs in the sectors (0, ±), (—1, ±) outside the 
curve of marginal stability. The corresponding trace of the monodromy matrix thus follows 


from (3.50) and is given by 


TrM^ = 


^2 cosh ^ -|- e {a+2ao)'^ 
(^2coshf + 


,(++-) 

(++-) 


if Im M > 0 (graph [^) , 
if Imu < 0 (graph [^). 


(3.51) 
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To compute the A-monodromy for the sector (H-), we can either cross four Stokes lines, 

or cross two Stokes lines and two branch cuts. The two computations are related by the 
relation (3.41). The former choice translates, in the case that the hrst Stokes line crossed is 
dominant, into the sequence of transition matrices 


with 




for [6b and 


iVl^2 


^2 —>■ 3 SrecSdom 

ao 

0 \ 

; e 

0 

) ’ 

e 

0 \ 

0 



:i-i c-i 


-^2-5-3 — 


N2^3 — 


ap+a 

e E 


CL j-^ — a 


0 

o.£)+a 

e 



(3.52) 

(3.53) 


(3.54) 


for 1^. The transition matrices upon shifting 0 hy tt are obtained by exchanging Sdom and 
Srec, which results in changing the sign of both a and a^) in the monodromy matrix. We 
thus obtain 


TrM^ = 


(^2cosh ifImM>0 

( 2cosh g±7« ) ifImM<0 

V " A+—) 


(graph 1^), 
(graph [^). 


(3.55) 


Upon inspection of £gure[^ one concludes that the results for Imu > 0 should be mapped 
to those for Imu < 0 via the map an ^ ao — this relation is in accord with the results of 
our computation. 


As was the case outside the curve of marginal stability, we again obtain multiple results for 
the trace of the monodromy matrix at a given fixed parameter u. The orientation dependence 
as reflected in the ± in the exponents in (3.51) and (3.55) is again resolved upon hxing the 
parameter pair {u, e): only one orientation is compatible with a given choice of e due to the 
relation (3.30). We will see in section 3.3.3 that the apparent difference between the results 

for the sectors (+ H—) and (H-) is accounted for by the jumping behavior of Voros 

multipliers. 


The 5-monodromy From the point of view of the Seiberg-Witten curve, a natural next 


task is to determine the i?-monodromy of the solutions to the differential equation (3.13) 


However, from the point of view of the g-plane on which the differential equation is for¬ 
mulated, the notion of H-monodromy is not meaningful, as the turning points qup/down are 
not identihed. The best we can do is ask how the wave functions obtained from Borel 
resummation of the WKB solution are related at qup/down- 


'4^ii.qdown^ t '0±(Q'iip) 


1 


's/ ^oddi^q Up) 
1 


exp ( ± 


rQup 


Sodd{q) dq 


\/ ^odd^q K 


:exp ( ± 


down ) 


^Qup 


'^Qdown 


+ 


)Sodd{q) dq 


'Qdo 


(3.56) 

(3.57) 
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In the final step, we have invoked the parity of the formal series Sodd- As both inside and 
outside the curve of marginal stability, an angle 6 exists at which we can connect qup and 
Qdown without crossing a Stokes line, we can elevate this relation to the level of the Borel 
resummed functions: 

(3.58) 


^:ki,Qdown} t ^ ^:izi,Qdowri} 5 


where the arrow indicates transport along a path connecting qdown to qup. On the Seiberg- 
Witten curve S/~, the quotient given by the identification q —)■ —q along the branch cut, 
equation (3.58) indeed describes the i?-monodromy of tki. 


3.3.3 Jumping phenomena for Voros multipliers 


It is possible to choose the normalization of WKB solutions such that their Borel transform 
remains locally invariant even as 6 traverses a critical value ig. This result implies that 
Voros multipliers must be sector dependent, and allows the determination of their transition 
behavior |T1] . Let 6 *o denote a critical angle at which a Stokes graph exhibits a single double 
Stokes line £q. Given a cycle 7 on the Riemann surface S (the double cover of S, see 
discussion below 3.32), the associated Voros multipliers (e^)_ and (e^)+ at 6*0 — 5 and 


60 + 6 , S a small positive constant, are related as follows M- 




(3.59) 


Here, 70 is a cycle on S whose projection onto S encircles io, with orientation chosen such 
that 

Ree*® i Xdz<0. (3.60) 

^70 

The intersection pairing (•, •) is chosen such that upon projection on E, the real and imaginary 
g-axis have intersection number one. 


Outside the curve of marginal stability Above, we labelled the sectors in between 
critical angles by an integer n determining the topology of the Stokes graph, and a sign 
determining orientation. The integer n increases with 6 *, as we will argue in appendix 
The double Stokes line l{n,±) which occurs at an angle 6 on the boundary between sectors 
{n — 1 , ±) and (n, ±) is homologous to the cycle ±(H + riH), and therefore the cycle 7(n,±) 
encircling it, with the orientation choice given by (3.60), is homologous to ^{2B + 2nA). 


By the property (3.59), the Voros multiplier associated to this cycle does not jump at the 


splitting of the Stokes line, i.e. 




(3.61) 


Noting finally that 


(7(n,±)5 T7 a) — 1 ) 


(3.62) 
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with 7 a of S a simple cover of the A-cycle on S, we can establish the independence of the 
trace of the monodromy matrix from the sector in which 6 lies as follows: 


TrMA 


-^{2na+2an)\ 


gf _j_ g“7 g^f («(l+2n)+2a£)) I 


/ (n,±) 


/ {n,±) 
-(7(n,±)’T7A) 

(n,±) 




(n,±) 




(n,±) 

-(7(n,±)’=F7A) 


(n-l,±) 


g±f g=Ff ^ g=F7(a(l+2(n-l))+2ao)^ 


/ (n-l,±) 


(3.63) 

(3.64) 

(3.65) 

(3.66) 

(3.67) 

(3.68) 


Inside the curve of marginal stability Let us discuss the transition from the sector 

(+ H—) to (H-) at Imw > 0 (figure to[^). The one at Imw < 0 will then follow upon 

the mapping a ^ an — a. The cycle 70 in this case is given by 70 = T2i?, where the — sign 
is for the case that the first line crossed in| 6 H is dominant. It follows that 


and hence 

TrMA 




„T- 


(+—) 




(3.69) 


± 


e' + e ‘ + e' 




) = (f±f(i + e±t”) + A?) 


_i_7m 
e £ 


1 + e 


=F- 


(++-) 


(++-) 


2a£)+Q _ 2a£)+Q g 

e £ + e £ + e+< 


(+—) 


(3.70) 

(3.71) 

(3.72) 


3.3.4 At the singular points k = ±1 


The e-neighborhoods of the singular points u = ±1 have received particular attention in the 
literature [ni[7iEniEiii2niiS2], as one is driven to these points at small real e and constant 


real characteristic exponent z/ (see the discussion around equation (3.10)). From the gauge 


theory perspective, m = ±1 are the points on moduli space where an extra state becomes 
massless and the effective gauge theory description breaks down. Nevertheless, we briefly 
touch upon this region in this subsection, and show how existing results for Tr Ma align with 
those we found above. 


The analysis presented in subsection |3.3.2| must be modified for the choice of modulus u = 
±1 -|- ef/, as the Stokes analysis depends on the function Qo, the leading term of Q{q,e) as 


introduced in (3.12), and is hence blind to the e U distance from the singular points. At these 
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points, the two turning points of Qo coalesce into a double turning point, and the connection 
matrices (3.39) and (3.40) are no longer valid. A calculation in the same spirit as the ones 
presented in subsection 3.3.2 was performed in [SO] at these points, proving a conjecture of 


PE! for the form of the characteristic exponents for this choice of u. The expression for 
TrM^ as cited e.g. in equation (2.14) of [29| for u = —1 can be expressed in terms of the 
periods a and ud by invoking the formulas (3.28), (3.8), and (3.33) of [29], yielding 


TrM^ = ^g-ia^g±i(a+ao)-^ 

f 2 cosh 

l 2 coshf+ 


(3.73) 

(3.74) 


with the periods a and evaluated ai u = —1 + eU G [—1,1]. This result coincides with 
the characteristic exponents for Imw > 0 inside the wall of marginal stability found above, 
for a particular orientation of the Stokes graphs. The two signs ± in the hrst line of this 
equation are due to choosing different branches for the term (— 

which emerges from equation (2.4) of [29]. A careful matching of conventions should allow 
reproducing the Imw < 0 results as well. 


4 Non-perturbative Effects in the 2d/4d Dictionary 


Our strategy in computing the instanton partition function exp of the gauge theory 

in PP was to compute the monodromy of the formal WKB solution to the null vector 
decoupling equation in the 62 —t 0 limit as a function of u, and to then determine m as a formal 
power series by imposing that this monodromy have no e-dependence. In this paper, using 
exact WKB methods, we have determined exponentially suppressed contributions ~ exp[— 
to the relation between u and the monodromy of the Floquet solutions to the differential 
equation. To analyze how these corrections manifest themselves in the instanton partition 
function, we wish to invert the relation between u and the monodromy, and then invoke the 
dehnition (2.25) of u in terms of T to determine exponentially suppressed corrections to the 
latter. 


The framework within which we will be performing these computations is that of transseries, 
a concept we will review in the next subsection before performing the calculation outlined 
above in subsection 14.21 


4.1 Transseries 

We will here consider the simplest class of transseries, which are formal power series in a 
hnite number of generators (see for instance USD - for our purposes, these consist e.g. of 
( 0 ) 

e, exp[—^], and exp[—^], organized term by term with regard to the obvious relation S> 
(with the relative size of the last two generators immaterial). To make sense of such series, 
we will require the formal power series in e that can be extracted at each exponential order 
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(0) ^\U J 

in exp[—and exp[—to be Borel summable. E.g., the transseries to the left of the 
arrow in the expression 

^ e-'^Sg{J2cmen ( 4 - 1 ) 

m m 

is to be associated to the analytic function in e (within a given sector) to its right. Transseries 
hence allow us, at the level of formal series, to distinguish between functions with the same 
asymptotic power series expansion, such as exp[—and 0. Furthermore, the relation S> 
lifts to the level of Borel transforms: an exponentially suppressed pre-factor cannot be 
compensated upon multiplication by the Borel transform of a formal power series, as a 
formal power series cannot capture the asymptotics ~ exp[—^]. 


When manipulating transseries, two questions naturally arise: is the result again a transseries, 
and can it be mapped to a function via Borel resummation (or a generalization thereof)? 
The former question is naturally easier to address than the latter, but also of limited useful¬ 
ness: subdominant contributions are only well-defined in the context of the map (4.1), and 
will generically depend sensitively on the integration direction 6, as we have seen. 


Transseries arise naturally in the context of exact WKB solutions. By removing the brackets 
(•), in (3.50), (3.51) and (3.55), we map our results for the trace of the monodromy matrix 
into transseries form. The coefficients of these transseries depend on u. The manipulation 
we wish to perform on these transseries is to solve them for u. As we will demonstrate in the 
next subsection, this is possible formally, and yields u as a transseries in the generators e, 
exp[—^], and exp[—A proof that this series is Borel summable is perhaps possible, 
combining information about the u-dependence of the coefficients (they are hypergeometric 
functions) and the growth behavior in e as follows on general grounds from WKB theory, 
but goes beyond the confines of this work. 


4.2 e-instanton corrections to the instanton partition function 


Since we are here interested in non-perturbative corrections to the instanton partition func¬ 
tion, we will consider the large u regime. Our goal is thus to solve the transseries equa- 
tioi]0 

eT + (4.2) 


underlying (3.50) for u expressed as a transseries in the generators e, e and e 


ao(nex) 


Recall that a{u, e) and aniu, e) are defined as formal power series in (3.23). In (4.2), we have 
expressed od as a function of a and e, which is always possible at large u. We will proceed 


by first solving (4.2) to express a as a function of and then plugging this relation into 


u{a, e) to obtain the desired result. With the ansatz 


d — (2ex ~\~ ^ 


m,n=l 


Cmrt.S 


- i (2m aex+2n au (aex,€)) 


(4.3) 


®We choose the sector (0, -I-) for convenience. The dependence of our result on the sector is a reflection 
of the shortcomings of transseries in the absence of a resummation result, as discussed above. 
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the first step amounts to solving consecutiveljj^ for the coefficients Cmn- Note that these 
coefficients will generically be formal power series in e, as they depend on derivatives of 
ao^a) evaluated at The hrst few terms are given by 


a = a. 


_ ^ g (2aex+2a/)(aex)) _|_ g ^ (4aex+2ai:>(aex)) _j_ —^3 ))c e ('^^ex+4a£)(aex)) _|_ ^ 

(4.4) 

The hrst terms of u{a, e) as a formal series in ^ and ^ can be computed based on the 
information provided in subsection 3.2 to be 


f(Osw, €) - ioL + ^ + ...+«" 




+ 


+ 


(4.5) 


where we have simplihed the expression at the expense of introducing yet another a-variable, 

\/2a 
2771 


Osw = The dual period a/? as a function of Osw is given by 


V2 . . 1 / , 1 15 

2^^®-d(®sw) f) — ^8®sw 2asw Sog-w H o 1“ xx ^ 1“ 


®sw 82 ag^ 


+e" 


3a; 


+ ^ + ... I +... 

sw ^sw 


(4.6) 


Plugging (4.4) into (4.5) yields the parameter u in terms of Oex, 

2aex T 2a£)(agx) 


u = 


(af^ - 2eae, exp 


dvr 2 


O (higher non-perturbative) 


) + 0(£=) 


(4.7) 


We can restore the scale A using 

e = ei/(y2A) and = a,tj{\f2K) (4.8) 


and hnd 


u = 


2eiarftexp 


2acft + 2a 


D,cft 


87r2A2' 

- O(higher non-perturbative) . 


ei 


(a.a/(y2A)) j 


) + 0(6?) 


(4.9) 


Integrating = —8Au with regard to A thus yields the following hrst ei non-perturbative 
correction to the instanton amplitude: 


= ^“cft log ^ + x-e? exp 


TT^ 


27r 


2 a eft + 2aD_eft(®cft/(v^A)) A 


+ ... . 


(4.10) 


Note that the expression (4.2) is invariant under the residual Z 2 gauge symmetry —>■ —a^ 


To proceed, we choose a particular sign for such that the exponential corrections in (4.3) 
are small. This choice corresponds to a gauge hxing of the residual symmetry. 


^One possible order is to solve for all m -I- t for consecutive i. 
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5 Conclusions and Open Problems 


We have seen that exact WKB methods strongly suggest that the general form of correction 
to the instanton partition function Zinst in an e-expansion be in terms of powers of exp[— 
and exp[—We computed such corrections to Zinst, the result arising upon formally 
solving a transseries equation for one of its parameters, u. It remains to be shown that 
the formal series thus obtained for the parameter u indeed corresponds to a transseries 
in the strong sense of subsection 4.1 We checked the consistency of our computation by 


demonstrating independence of our result from the integration direction 9 of the inverse 
Laplace transform. The specihcs of the transseries form do depend on this parameter. 


The techniques used in this paper can equally well be used to compute non-perturbative 
corrections to superconformal gauge theories, such as A/" = 2*. The modularity constraints 
in such theories may allow for the extraction of more detailed non-perturbative informa¬ 
tion. 


More broadly, we relate Z^st to the null vector decoupling equation by taking the limit 62 —)■ 0 
in parameter space for which no worldsheet description of the topological string theory is 
yet available. Filling this gap, or adapting our methods to apply away from this limit, will 
be an important step towards the further physical interpretation of our results. 
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A Local Analysis of Stokes Lines 


Consider the part of parameter space where |m| S> 1. In this case, the two turning points of 


Qo in equation (3.14) lie far from the real axis, 

Qtp ~ log 2u + 27m , 


(A.l) 


and the cosine can be well approximated by an exponential. For a turning point in the upper 
half-plane. 


Xdq ^ j \j-e-X<itp+x) 


(A. 2 ) 


The angles (px at which the Stokes lines emerge from the turning point are determined to 
hrst order by the condition 


1 2 1 

-71 + (px + 4>u = ^{-9 + rni - -71 + (pu), neZ, 


(A.3) 
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with u = 


lines hence emanate in the directions 


= -|0+|7r(n + |) 


Note that increasing 6 results in (j)x decreasing, the Stokes lines hence rotating anti-clockwise. 
This implies that when decreasing an angle from a critical value, the two simple Stokes lines 
into which the double Stokes line splits swerve to the left, as seen from the turning point from 
which they emerge, whereas when increasing it, they swerve to the right. This is a property 
we used in subsection |3.3.3[ Let us now consider u > 0 and 6^ = 0 for illustration purposes. 


In this case, the Stokes lines occur at angles (px = ^tt, |7r. Which n corresponds to 

each of these three angles depends on the choice of branch cut, i.e. the choice of interval 
for the phase of the argument of the square root. Note that the pattern of Stokes lines is 
hence independent of the choice of branch cut, whereas the orientation of the Stokes lines 
does depend on this choice: for n even, the Stokes line is dominant, for n odd, it is recessive. 
To have e.g. the Stokes line in the negative imaginary direction have different parity from 
the other two, we can choose the argument of the radicand to lie in the interval [0,27r], 
corresponding to a branch cut in the positive imaginary direction. The values of — ^vr -1- (px 
(the phase of the argument of the square root) that fall into this interval are ^tt, tt, |7r, with 
TT corresponding to the negative imaginary direction. The corresponding values of n are 
1,2,3, hence the distinguished Stokes line is dominant]^ For a turning point in the lower 
half-plane, we obtain 


/ 


A dq 


r\j 



i(qtp+x) _ ^ 


dq 


rsj 


2i 

3u 


{ixu) 


3 

2 

1 


(A.4) 


hence with the same choice of branch cut, (px = corresponding to n = 0,1,2. 

The Stokes line pointing in the positive imaginary direction is hence recessive. 


B Numerical Results 


In this appendix, we will compare results for the Mathieu characteristic exponent obtained 
numerically with those attainable by evaluation of the formulae determined in subsection 

ESI 

B.l Numerical WKB perturbation theory 

In a hrst numerical experiment, let us illustrate the efficacy of WKB perturbation theory. 
Consider the following values of u and e: 

e = ^, u = 2 + 2i. (B.l) 

®Note that there is a small subtlety here due to the nature of branch cuts. The branch of the square 
root is specified by specifying the interval within which the phase of its argument lies. If the argument is a 
function, we usually specify this interval for the variable of that function. But even in the simplest case of 
a linear function az, this requires specifying the phase of the factor a exactly, i.e. not only up to multiples 
of 2tt. This is why we specify the branch cut not in terms of <px, but in terms of the argument of —ix. 
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The value of e is small compared to 1 and the value of u lies outside the curve of marginal 
stability. The values of a and 2a d — a si leading order in perturbation theory ar^ 

aW/e ~ 41 - 97i , {2aD - ~ -90 + 19i. (B.2) 


Non-perturbative effects are hence heavily suppressed, and the Mathieu characteristic expo¬ 
nent should be well approximated by the a-period of the differential A and its higher order 
corrections in e. Indeed, the Mathematica routine for the numerical approximation to the 
Mathieu exponent gives the value (up to 20 digitsj^ 

p = 41.209556582920410400 - 97.333703725326654896Z. (B.3) 


The successive approximations using the WKB formulae (3.26) in e-perturbation theory 
are 


= 41.20893 - 97.33334i, 
a(2)/e = 41.20955667 - 97.33370309*, 

= 41.2095565837- 97.33370372575*, 

= 41.2095565829175 -97.3337037253288*, 
a(8)/e = 41.209556582920406 - 97.33370372532663*. (B.4) 

Each iteration enhances the accuracy by approximately two digits, which is what one would 
naively expect of a perturbation series in e^ = 10“^. These numerical results demonstrate 
the usefulness of WKB perturbation theory, even in its asymptotic form. 


B.2 Numerics and the Stokes corrected formula 


In this subsection, we will consider examples for which the non-perturbative effects computed 
in subsection 3^ become numerically signihcant. The evaluation of the formulae derived 
there will be sector dependent, as we will approximate Voros multipliers by the hrst terms 
in their asymptotic expansion. We will work inside the curve of marginal stability. Of the 
two sectors that occur here, we will observe that one yields the better approximation to 
the characteristic exponent, indicating that the Voros multiplier in this particular sector is 
better approximated by its asymptotic series. 


Let us consider the values 


1 

^ “ 5 


u = -e *10 


(B.6) 


We shall compare the value for the trace of the monodromy matrix determined by Mathe¬ 
matica [53] to those which follow from the numerical evaluation of the perturbative WKB 


®In this appendix, the symbol for a formal power series with superscript will indicate that power series 
numerically evaluated up to the n-th order term. 

lOWe record the Mathematica value times —in to render normalizations uniform. 

^^Here and in the following, we do not round the numerical values. 
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result and the exact WKB result in the sectors (+ H—) and (H-), to 0*^ and 4^^ order 

in perturbation theory. 


y^^Matheniatica 

1 

2 cosh -j 


2 cosh —h e' 
e 

2aD — a 


2 cosh 


+ e* 


( 0 ) 


1 

2 cosh -j 


(2 cosh® +e^(2az,-a)j 


(4) 


3078.40577 - 11972.57629*, 
3313.34282 - 11850.61433*, 

2810.06097 - 12265.24243*, 
2810.05976 - 12265.24153*, 

3582.79033 - 11581.22299*, 
3078.40503 - 11972.57314i, 


2 cosh 


2aD — Cl , 1 , 

-h 


(4) 


= 3078.40377 - 11972.57226*. 


(B,6) 


The formulae for the two sectors (+ H—) and (H-) yield numerically close results, as we 

are approximating 


2 cosh —h e' 
e 

o — a 

2 cosh- 


{2au-a) 

) 

~ + e 

/ 

+++-) 


■ + 

) 

~ + e 

J 

'(+-) 



+ e“ 


(B.7) 

(B.8) 


with 

a(4)/e = 9.40284 + 17.5788*, (2aD - a)(^Ve = 6-45897 - 15.0481*. (B.9) 


To be able to better discriminate numerically between the two equations (B.7) and (B. 


we shall next choose a value of u at which |2a£) — a| is small, and then choose the phase of 
e to eliminate the real part of a. 


To this end, let us choose 


Then 



(B.IO) 


oW = 4.23320+4.3172; = 6.04636e‘’'l“-^“'^>, (2oD-o)<'‘Ve = 0-04966-0.28854i. (B.ll) 

Choosing e to rotate away the phase of a. 


e = 



(B.12) 
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we obtain 


TrM^ 


Mathematica 


9.84414, 

-1.46994, 



10.03070, 


10.85260 + 0.67809i, 


-1.46174, 


(2 cosh 


9.84414, 



10.66350 + 0.68251i. 


(B.13) 


For this example, the approximate evaluation of our exact WKB formulae in the sectors 

(H-) and (+ H—) yields appreciably different values, and one sector, (+ H—), yields the 

better approximation. This indicates that the Voros multipliers are better approximated by 
their asymptotic expansion in this sector. 

C Periodicity, Determinant, and Numerics 

In this section, we review the determinant formula for the exact periodicity, an efficient 
way to evaluate the determinant, and how the resulting numerics may improve on a built-in 
Mathematica [53] evaluation of the Mathieu characteristic exponent. 


C.l Hill’s method 


Hill’s method yields an exact formula for the Mathieu characteristic exponent v in terms of 
the parameters in the Mathieu equation (see e.g. [5l| and references therein). 

The derivation of this formula proceeds as follows. The power series 



(C.l) 


provides a formal solution to Mathieu’s equation (3.13) if its coefficients satisfy 


2C2r-2 + (6^(2r -|- z/)^ — 4M)c2r- + 2C2r+2 — 0 . 


(C.2) 


By dividing all coefficients in this recursion relation by the coefficient of the second term, the 
determinant A(z/) of the matrix underlying this linear set of equations for the coefficients Cr 
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becomes convergent 0. For a non-trivial solution to this infinite set of equations to exist, 
A(z/) must vanish. The determinant is invariant under n —)■ 2n ± z/ (with n integer), and is 
therefore even and periodic in u with period 2. It has simple poles at (2r + z/)^ — 4tze“^ = 0, 
and it tends to 1 as z/ —>■ oo. This behavior determines, by Liouville’s theorem, the form of 
A(z/) up to a constant [M] 


1 - A(z/) = 


sm“ ^ — sin^(7rA/ue^) 


Imposing A(z/) = 0 then gives rise to the following constraint equation on z/: 


— = A(0) sin^ TiVue ^, 


sm 


or equivalently 


cos TTz/ = 1 — 2A(0) sin^ vrV ue~‘^ . 


(C.3) 


(C.4) 

(C.5) 


C.2 Linear Recursion 

The determinant A(0) can be computed efficiently as it satisfies 

A(0) = detv4, (C.6) 

where A is a tri-diagonal matrix |H1]. We denote the determinant of the submatrix of size 
{2i -f 1) X (2z -|- 1) acting on (c_ 2 i,..., c_ 2 , Cq, C 2 , ..., C 2 i)'^ as Aj. In the limit z —)■ oo, we 
recover A. The determinants Aj satisfy a linear recursion relation (proved by computing 
minors) 


Aj (1 Q;j)Aj_]^ o^z(l Q;j)Aj_2 T Q;jCrj_^Aj_3, 


(C.7) 


with 


and initial terms 


1 

A{n‘^e^ — u){{n — 1)^6^ — u) ’ 


Aq 

Ai 


A 2 


1 , 

(1 

det ^0 

V 0 

/ 1 


det 


6 

0 

0 


V 0 


6 0 \ 

1 ^0 , 
e-2 1 / 

^40 0 0 \ 

16 0 0 

6 16 0 

0 62 1 62 

0 0^-4!/ 


(C.8) 


(C.9) 


39 










where 


2 

_ 4 ^ 


(C.IO) 


= 


Our implementation of this linear recursion in Mathematica demands a computation time 
that grows linearly in the size of the matrix A. The reduced cost is due to the fact that A 
is a tridiagonal matrix. 


C.3 The Determinant, Mathematica, and WKB 


Mathematica provides a convenient numerical algorithm for hnding the Mathieu character¬ 
istic exponent for any complex parameters of the Mathieu equation. However, the algorithm 
is a black box. In our numerical experiments, we have found that the algorithm must be 
used with caution, in particular at small values of |e|, as the following example demon¬ 
strates. We will compute the characteristic exponent in a particular case using three meth¬ 
ods, namely Mathematica, WKB perturbation theory, and the numerical algorithm based 
on Hill’s method describe in subsections 1C.II and 1C.21 


We consider the values 

e = 10-\ M = 6e*t, (C.ll) 

and have Mathematica compute the parameters 

aM = 4u/e^, qM = 2/e^ (C.12) 

and the characteristic exponent 

N[—i7rMathieuCharacteristicExponent[aM, Q'm], 20]], (C.13) 

yielding the result 

301.3467972919577793 - 1507.3089399330755493z. (C.14) 

The WKB approximation to fourth order for the a-period yields 

301.7573905535926129 - 1507.2622454539553016z. (C.15) 


At small e, we expect to obtain a better approximation to the characteristic exponent than 
the comparison of (C.14) and (C.15) would suggest, given that a/e ~ 302 — 1507z and 
2aD — a —2020 -|- 770i imply highly suppressed non-perturbative corrections. 


To check the Mathematica result, we programmed a numerical algorithm based on Hill’s 
method. For a matrix of size 2 x 10®, we found the result 

301.75739055567594974 + 0.70222826914545431z. (C.16) 

This agrees with WKB perturbation theory to high order ~ the difference in the imaginary 
part is an integer multiple of 27r, which corresponds to the ambiguity in the Mathieu charac¬ 
teristic exponent. Thus, for small e, the perturbative WKB approximation is more reliable 
than the (unspecified) Mathematica algorithm. It agrees with an alternative numerical al¬ 
gorithm which has linear cost in the size of the (sparse) matrix used in the evaluation of the 
determinant. 
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